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We study the three-dimensional bosonic t-J model, i.e., the t-J model of "bosonic electrons" 
at finite temperatures. This model describes a system of cold bosonic atoms with two species in 
an optical lattice. The model is derived from the Hubbard model for very large on-site repulsive 
interaction between bosons of same species (hard-core nature) and also strong correlations between 
different species. The operator B^cy for an atom at the site x with a two-component (pseudo-) 
spin cr(= 1,2) is treated as a hard-core boson operator, and represented by a composite of two 
slave particles; a spinon described by a CP^ field (Schwinger boson) z^^ and a holon described 



by a hard-core-boson field 



as iJxa — 



At. 



is then expressed by a pseudo-spin, which 

We then have a 



is, in turn, represented by another CP^ (pseudo) spinon Wj,,, as (j)x ~ wl.2Wxi- 
double-CP'^ representation of the model by Zxa and lOa;,,. By means of Monte Carlo simulations of 
this bosonic t-J model, we study its phase structure and the possible phenomena like appearance of 
antiferromagnetic long-range order, Bose-Einstein condensation, phase separation, etc. They should 
be compared with the possible experimental results of a recently studied boson-boson mixture like 
*^Rb and *^K in an optical lattice. 

PACS numbers: 67.85.-d, 11. 15. Ha 



I. INTRODUCTION 

Cold atoms [i, 0] trapped in optical lattices are one 
of the most interesting systems in condensed matter 
physics. They are flexible because one can experimen- 
tally adjust their characteristics such as statistics and 
density of atoms, dimensionality of the system, signature 
and strength of interactions, etc. 

Not only systems of single kind of atoms, but also sys- 
tems of atoms of two species like a mixture of bosons 
and fermions are studied experimentally [3|. Recently, 
a boson-boson system of *^Rb and ^^K in a three- 
dimensional (3D) optical lattice has been studied to 
produce Bose-Einstein condensate [J]. The interactions 
among these atoms are also controllable Q. 

A standard model of cold bosonic atoms with two 
species may be the bosonic Hubbard model. The two 
species of bosons may be described by a (pseudo-)spin 
s = 1/2 degrees of freedom. For very large strong re- 
pulsive interactions between atoms of same species, one 
may treat each boson as a hard-core boson (HCB). Thus, 
usual electrons of Hubbard model are to be replaced here 
by HCB with spins. 

From this HCB Hubbard model, one can derive the 
bosonic t-J model as its low-energy effective model for a 
large on-site repulsion between the opposite spins (differ- 
ent species) and small hole concentrations p,0|. 

By using this bosonic t-J model, Boninsegni and 
Prokof'ev[8] studied the interplay of magnetic order- 
ing of pseudo-spins and Bose-Einstein condensation 
(BEC)/superfluidity (SF) of bosonic atoms. By quantum 
Monte Carlo (MC) simulations, they studied the low- 
temperature (T) phase diagram of the two-dimensional 



(2D) model for the case of anisotropic spin coupling 
Jx.y — ctJz, a < 1 and J^ = J < t, and found the coexis- 
tence region of antiferromagnetic (AF) order and SF as a 
result of the phase separation (PS) of hole-free (AF) and 
hole-rich (SF) phases. 

In the previous paper 9] we studied the bosonic t-J 
model with the isotropic coupling {a = 1) in the slave- 
particle representation of operators for atoms. The use- 
fulness of the slave-particle representation in various as- 
pects has been pointed out for the original fermionic t-J 
model[l3| . We expect that similar advantage of the slave- 
particle picture holds also in the bosonic t-J model[ll|. 

In the slave-particle representation, the bosonic opera- 
tor Bxa- for atom at the site x and spin a = 1, 2 is viewed 
as a composite of a spinon z^a and a holon (jj^ [9] , 

B^a = (l^lzxa-, (1-1) 

where Zxa- is the CP^ spin field (Schwinger boson) and 
(j)^ is the HCB. 

In Ref . 0] , we replaced these HCB operators of holons 
(j)x by the Higgs field with a definite amplitude, 



(t)x -> \/Se:xip{iipx), 



(1.2) 



where S is the average density of holes (holons) S = 
{4>l4'x) and (fix is the phase degrees of freedom. This 
replacement is an approximation to ignore the fluctua- 
tion of amplitude of holons assuming their homogeneous 
distribution. By the MC simulations we have obtained 
a phase diagram and various correlation functions of the 
3D model at finite temperatures (T > 0) and the 2D 
model at zero temperature (T = 0). In both cases, we 
found the coexistence region of AF order and SF. 

It is then interesting to relax the above assumption 
of homogeneous distribution of holes and consider the 



possibility of PS. With an isotropic AF coupling, it is 
reported [l2| that the ground state of the 2D bosonic t-J 
model is spatially homogeneous without PS for J/t < 1.5. 
In this paper, we shall study the 3D bosonic t-J model 
at finite T in the slave-particle representation without 
the above mentioned approximation. That is, we treat 
the holon variables (p^ as genuine HCB instead of the 
compact Higgs field exjp(iipx). We express (px as another 
CP^ (pseudo) spin field Wxtj {rj = 1,2) via pseudo-spin 
SU(2) operator as 

4'x = wl^Wxi- (1-3) 

We study this double-CP-^ system of Zxu and Wxt^ by MC 
simulations. We examine its phase structure, various cor- 
relation functions, and possible PS, and so on. 

The present paper is organized as follows. In Sec. 2, we 
introduce the double-CP^ representation of the bosonic 
t-J model in the 3D lattice at finite T's. We define two 
versions of the model. Model I directly derived from the 
t-J model, and Model II, a simplified version of Model I. 
They have different weights in spin stiffness. By compar- 
ing the results of these two models, one may obtain fur- 
ther understanding of the interplay of holons and spinons. 
In Sec. 3, we exhibit the results of the numerical study and 
the phase diagram of the simplified mode. Model II first. 
We calculated the specific heat, the spin and atomic cor- 
relation functions. From these results, we conclude that 
there exists a coexisting phase of AF long-range order 
and SF in a region of low-T and intermediate hole con- 
centrations. In Sec. 4, we consider Model I and present 
the MC results. We compare the phase structures of 
the two models. We also study the PS of both models. 
Section 5 is devoted for conclusion. In Appendix A, we 
derive the expression of the HCB operator (px in terms of 
another CP^ operator Wxr)- 



II. MODEL I AND MODEL II 

Let us start with the bosonic Hubbard model of HCB's. 
Its Hamiltonian is given by 



i?i 



Hub 



-* E i^l 



X + fJ,,(7-^X(T ■ 



H.c. 



U^nxinx2, 



"'Xa ^-'rrr-'-'X 



(2.1) 

where Bxa- is the HCB operator [Tj, [Ijl to describe the 
bosonic atom at the site x of the 3D cubic lattice and 
the spin a[ = l(t),2(;)]. ^(= 1,2,3) is the 3D direction 
index and also denotes the unit vector. Bxa satisfies the 
following mixed commutation relations for HCB's[15|. 



Bxa,B^^, 



' t Y-t^xa-i J^xa' 



0, 



Bxrr,B, 



ya' 



= [Bxa,Bya,]=OioTx^y. (2.2) 
We note that ifnub has a global SU(2) symmetry under 

Bxa -^ B'x, =Y,9aa-Bxa-. 9 £ SU(2). (2.3) 



Actually, the first term is manifestly invariant. To show 
that the second term is also invariant, we introduce the 
atomic number operator. 



nxi +nx2, 



(2.4) 



which is also invariant under (12.31) . Then the second term 
is expressed sorely by this invariant quantity fix as 



nxinx2 = linxi 



1 



-nx2 



(nil + ^D 



-(nl -fix), 



where we used the relation fiz 



(2.5) 



which holds be- 



cause the eigenvalue of fixa is or 1 due to HCB. 

From this Hubbard model one may derive the bosonic 
t-J model as its effective model for strong correlations, 
i.e., for large U, and small hole concentrations. The 
Hamiltonian H of the 3D bosonic t-J model is given by 



H = -iJ2 (^1+M,A.+H.^ 



+ <-' / ^ I ^X-\-fl ' ^X 
X,fi 



''X-\-fJ,"'X I ; 



Bxa^il-flxa) Bxa, (1 = 2,2=1) 



Sx = —B^aBx, {S : Pauli matrices). (2.6) 



The derivation of (|2.6p is achieved just by following the 
steps developed in the theory of high-Tc superconductiv- 
ity to derive the fermionic t-J model from the standard 
Hubbard model. Actually, the second-order perturbation 
theory in terms of small parameter t/U produces (j2.6p 
with the relation J = 4t^/[/[l^. The role of Pauli prin- 
ciple for fermions is played here by the hard-core nature 
of bosons. 

Because of the hard-core nature of Bxa, {BI.^Y = 0, 
the states of double occupancy of the same-spin bosons 
at each site are excluded. Furthermore, as we consider 
the case of strong repulsion between two bosons with 
opposite spins at the same site, we impose another con- 
straint that excludes from the physical space the double- 
occupancy state of two bosons with the opposite spins at 
each site. This implies the states -BL-BL |0) are excluded 

from the physical space[l3|. The tilde operator Bxa re- 
fiects this fact. These two constraints are expressed for 
the normalized physical states |phys} as 

(phys|n:r|phys) = (phys| y^ gt^^^rcrlphys) < 1. (2.7) 



The Hamiltonian (|2.6|) respects SU(2) symmetry 
To see it, it is helpful to rewrite Bxa = (1 ~ nxa)Bxa = 
(1 - nx)Bxa due to Bl^ = 0. 

As mentioned, the fermionic counterpart of the present 
bosonic model, the (fermionic) t-J model, is a canonical 
model for strongly-correlated electron systems like the 



high-temperature superconductors. The 2D fermionic t- 
J model has been studied intensively by means of vari- 
ous methods since the discovery of the high-T supercon- 
ductors, although its phase structure etc. are still not 
clarificd[18]. The slave-particle approach like the slave- 
boson or slave-fermion representations has provided us 
with an intuitive way of description of the t — J model 
in the mean field theory, in the charge-spin separated 
state[l^, and so on. It is natural to expect that the 
slave-particle approach is useful also in the bosonic t — J 
model. 

In the slave-particle representation, B^a is expressed 
as 



_D-r 



(2.8) 



(px represents the annihilation operator of the hard-core 
holon, and satisfies the HCB algebra of single species, 



''xtYx] 



/>tl - 



+ 



1, 



J+ 



-0, 



f>l]- 



bx,<t>y]^QioTx^y. (2.9) 

o-xa represents annihilation operator of the bosonic 
spinon carrying s = 1/2 spin. Their commutation re- 



lations are 



,t 



^xy^aa' 



[axcr,aya'] = 0. (2.10) 



The physical-state condition (|2.7I) is replaced by the fol- 
lowing equation in terms of the slave-particle operators 
axa- and (h^, 



Cyi al„o,xa + (t)i(l)x)\phys) = |phys}. 



(2.11) 



Meaning of Eg. (|2. lip is obvious. There are three physical 
states for each x; one-holon state (holon number (/)J,0x = 
1), corresponding to the state with no atoms, and the two 
one-spinon states (one with spinon number aliUxi ~ 1 
and the other with ax20.x2 = 1) corresponding to the 
two one-atom states. The correspondence among these 
physical states are given as follows: 

Bxa\0)=0, 

(j)x\veic) = a^ECT I vac) = 0, 

|0) = 0![,|vac), No — atom state 
BI^\0) — a^g.|vac). Two one — atom states. (2.12) 

In terms of the slave-particle operators, the Hamilto- 
nian (12.61) becomes 



H = -t 



X,:Lfl.(T 



(i^ui 



:tf-i,(T 



^xa^x 



±^l 



4^ 



[{a^aa)x+f, ■ {a^ffa)x - ia^a)x+^^ia^a)x] , 



(2.13) 



where we write (a^a) 



E<TaLax<T, {a''aa). 



The constraint (I2.11[) is solved as 

axa = (1 - 4>l<l)x)Zxa, (2.14) 

where Zxa is the CP"'^ spin operator that satisfies the usual 
commutation relations for bosons and the following CP"'^ 
constraint. 



/ , ZxfjZxa — 1- 



(2.15) 



To derive Eq. p.l4|) we have used the identity 

(l-40,)2 = (l-40,), (2.16) 

because the eigenvalues of 0J.0X are and 1. 

In Appendix A, we show that the HCB operator (^^ 
can be exactly expressed in terms of another CP^ op- 
erator Wxr) (j] = 1,2), which satisfy the usual bosonic 
commutation relations and the constraint. 






(2.17) 



The expression of 4>x is then given by 

<l>x^w\iWx\- (2-18) 

Then we employ the path-integral expression of the par- 
tition function 



Z = Trexp(-^iJ), 



(2.19) 



at T with /3 = ll{kBT). For this purpose we introduce 
two sets of CP-^ variables, Za;cr (t) , ui^^^ (r) at each site x 
and imaginary time r G [0,/3 = {k-QT)~^]. (Hereafter we 
often set the Boltzmann constant fcs to unity.) They are 
complex numbers satisfying 



'^Zxc{t)Zxo{t) = 1, 



(T=l 

2 



'^Wxn[T)Wxr){T) = 1. 
^=1 



(2.20) 



Then the partition function is given by means of path 
integrals over Zx(t{t) and Wxri{T)- 

To proceed further, we make one simplification by con- 
sidering finite-T region, such that the r-dependence of 
the variables in the path integral can be ignored keeping 
only the zero modes, i.e., Zx(t{t) -^ Zxa, etc.[9|. To study 
the phase structure at finite T is very important because 
it summarizes the essential properties of the system. Be- 
sides it, finite-T phase diagram gives a very useful insight 
into the phase structure at T = 0, i.e., if some ordered 
states are found at finite T, we can naturally expect that 
they persist down to T = 0. 

Let us present the path-integral expression Zi of the 
partition function (|2.19p of the 3D model at finite T's. 



Z\= j W [dzxdwx Y\ dUxf^] exp(yli - /i^ ^ (f'x^'x 



(2.21) 



where the suffix I has been attached because we call this 
model Model I (We shall introduce Model II later). The 
action Ai on the 3D lattice is given by 



Ah = 



C3 



.x,^,a 



Px = l~ (j)x(t)x, 
Zxl = Zx2, Zx2 = ~Zxl {Zx = io-2^*), 

4>x = Wx2Wxi, Uxii = exp(i6'a;p), (2.22) 

and the integration measure is 

//"OO /"OO 

dzx = / dzxi / dzx2 S(^ZxcrZxa- " l) etc, 
J —00 J — 00 „ 



dUxfj. 



— 00 
27r 



dO^fj, 

~2^ 



(2.23) 



Here we have used the same letter (j)j; for the complex 
variable and the operator in Eq. p.Qp because no confu- 
sions arise. According to Ref. 9] , we have introduced the 
U(l) gauge field Uxfi = exp(i02;^) on the link {x,x + fi) 
as an auxiliary field to make the action in a simpler form 
and the U(l) gauge invariance manifest. 

The term —fJ-c'^x^^^^ with the (minus of) chemical 
potential ^c has been introduced to control the hole den- 
sity p, 

p{ci,C3,fJ.c) = j^^ifjyx'px), (2.24) 

X 

where N — J^x ^ ^^ ^^e total number of the sites. As 
indicated, p{ci,C3, pc) is a function of ci,C3 and p-c- We 
are interested in the case that the hole density p takes 
a constant value S as ci,C3,pc sue varied, because each 
material has a constant hole density, S. Then, to obtain 
a physical quantity for given ci,C3 and 6, we determine 
Pc so that the relation 



p{ci,C3,pc) = 6, 



(2.25) 



holds. It implies to determine /Zc(ci,C3,(5) as a function 

of Ci,C3,5. 

Because the equation 



1 



(2.26) 



holds for the physical states, there is a relation among 
the atomic density n and the hole density p as 



-(^n,) = l-p. 



(2.27) 



as expected. 



The action Ai of ((2?22l) has the global SU(2) spin sym- 
metry of (1231), 

Zx -^ gzx, 9 e SU(2) 



'^X-\-fJ,^X 



^x-\-}.i99^x Zx-\-fiZxi 



^x-i-fi^x ~^ clet^ • Zx-\-^Zx — Zx-^i_iZx- yZ.Zo) 

We note that the present model may describe not only 
s = 1/2 case but also the case with an arbitrary s. In 
this case, Ea. (|2.20p becomes as[2Q| 






25, 



(2.29) 



but the above normalization factor 2s can be easily ab- 
sorbed into the parameters ci and C3 in the action. 

The action Ai is also invariant under a local (a;- 
dependent) U(l) gauge transformation. 



Zx(7 ^ G Zxai 



„i^x 



Uxf, -^ e 



+"C/.Me 



-lA^ 



Wxl 



e 2 Wxl, Wx2 -^ e' 



Wx2, 



(2.30) 



where Xx is an arbitrary function. The complex vari- 
able Bxa 21| for the gauge-invariant bosonic atoms is ex- 
pressed in terms of z^ and w,, as 



Bxa = (fixZxa- = WxlWx2Zxa 



(2.31) 



The parameters ci and C3 are related with those in the 
original t-J model as[9[ 



J/3 for d » 1, 

"^ ^ ^ (2J/3)i/2 for ci « 1, 

C3 ~ t/3. 



(2.32) 



Note that C3 here has no extra factor S compared with 
C3 defined in Ref. [9]. 

This model (|2.22p . which we call Model I, is examined 
in Sect. 4 as announced. Before that, we study a simpli- 
fied model, Model II, which is defined by setting P^ — 1 
in (|2.22p . The partition function Zu of Model II is then 
given by 

Zu^ /IT [^^^^^^ n '^^^m] exp(Aii - Pc'^<[>x4'x), 

X fl X 

Aii^ A^ + Ai,, 



fiZxa ~r C.C 



(2.33) 



The main reason to study Model II is to clarify the effect 
of the projection operator to the hole- free states, Px, by 
comparing the results of the two models. We expect that 
the AF ordered state appears with stronger signals in 
Model II than in Model I because the assignment Px = 1 
in Model II lets the AF coupling between the nearest- 
neighbor (NN) spin pairs at (a;, x + p) survive even if a 
hole occupies the sites x and/or x + p. Here we note 
that, in the path- integral formulation, the variable Zx(t) 



is defined for all (x, r) even if the site x is totally occu- 
pied by a liole[22|- Due to the short-range AF configura- 
tion, such Zx reflects a nearby spin orientation. In other 
words, Model II effectively describes a doped AF Heisen- 
berg model with, e.g., next-NN exchange couplings that 
enhance AF long-range order. 



III. RESULTS OF MC SIMULATIONS OF 
MODEL II 

In this section, we present the results of MC simula- 
tions for Model II of (|2.33|) . For MC simulations, we con- 
sider a 3D cubic lattice of the size N = L^ (L up to 36) 
and imposed the periodic boundary condition. We used 
the standard Metropolis algorithm with local updates. 
Average number of sweeps was 12 x 10"', and average 
acceptance ratio was about 40%'^70%. 



A. Phase structure 

Let us first discuss the phase structure. To this end, 
we measured the internal energy [/, the specific heat C 
and the hole density p defined as 

1 



iV^^)' 



u 



(3.1) 



In Figdlwe show the obtained phase diagram in the C3- 
ci plane for various values of p = 5. The phase transition 
lines in Fig [T] were determined by the measurement of U 
and C in (|3.1I) . We found that there are four phases listed 
in Table 1. The nature of each phase was confirmed by 
measuring various correlation functions, which we shall 
discuss in the following subsection in detail. 

Let us present a qualitative understanding of the phase 
diagram Fig IT] For this purpose, it is convenient to in- 
troduce an 0(3) spin vector tx made of spinon z^. 



ZxCTZx^ ^x ' ^x — ^j 



(3.2) 



where summations over spin indices are understood. The 
short-range (SR) spin correlation £x+ij,-^x is expressed by 
the SR CP^ amplitude (such as Zx+fj.Zx) as 



^a:+/i ■ ^x 



2|£ 



x-'rfi'^x I 



-2|£ 



a:+/i^x I 



1, (3.3) 



AF: Antiferromagnetic phase 

PM: Paramagnetic phase 

AF-I-SF: Phase of antiferromagnetism and SF 

FM-I-SF: Phase of ferromagnetism and SF 

Table 1. Various phases in Fig[T] 
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FIG. 1. Phase structure of Model II in the C3-C1 plane for 
various 5. The meaning of abbreviations for each phase is 
listed in Table 1. The phase transition into the FM+SF phase 
for i = 0.05 and the phase transition between the AF and 
AF+SF phases for 5 = 0.15 (solid curves) are of first order. 
All the other transitions (dashed curves) are of second order. 



where Zx satisfies 



1- |2 I I ~ |2 _ 1 

\^x+fi^x\ "r \Zx+^Zx\ — 1. 



(3.4) 



So the SRAF configuration (.x+fi ■ ix — —1 corresponds 
to \zx+iiZx\ — 1, and the SRFM configuration ix+ii ■ ^x — 
1 corresponds to \zx+fiZx\ — 1. The c\ term A'^ of 
(|2.33|) controls the spin fluctuations and the AF order 
(|za;+^Z3;| ~ 1) is generated for sufficiently large values of 
ci. The C3 term Ay, of ^?F^, (I05|) controls the rate of 
hopping of holons, (j)x+^(f>x, accompanied with a SRFM 
spinon amplitude Zx+fj.Zx- The SF of holons is generated 
for sufficiently large values of C3 if there exists a sufficient 
amount of FM spin amplitude {zx+fj,Zx)- Therefore, for 
sufficiently large values of C3, we expect a coherent SF 
of holons accompanied with a (SR) FM order. As the 
results in Fig[T] show, there exists the coexisting phase 
of AF and SF for intermediate and large hole density. 
More details of the spin configuration in this phase will 
be given in the following sections. 

In Fig 12] we show the phase diagram drawn in the 
S — 1/ci plane for several values of C3/C1. The value 
for vertical axis, 1/ci, and the ratio C3/C1 are expressed 
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FIG. 2. Phase structure of Model II in the S-l/ci plane 
(1/ci - ksT/J) for C3/ci(- t/J) = 2.0,3.0,4.0. The abbre- 
viations for each phase are listed in Table 1. The phase tran- 
sition between the AF and AF-I-SF phases for ca/ci = 3.0, 4.0 
drawn in solid curves are of first order. All the other transi- 
tions drawn in dashed curves are of second order. 



for ci >> 1 from (|2.32p as 



fcf 



t 
1' 



(3.5) 



ci J/3 J ' ci 

So Figl5] may be interpreted as the phase diagram in 
the 5 - T planedl]. The Neel temperature T^{5) (the 
critical temperature of the AF transition) decreases as 
holons are doped (as 5 increases). The rate of reduc- 
tion —dT^{S)/d5 increases as t/J increases because the 
movement of holons becomes significant to destroy AF 
order. The critical temperature Tsf{6) for SF of holons 
increases as 6 increases as expected. The rate of increase 
dTsF{S)/d6 also rises as t/J increases, because the move- 
ment of holons favors their homogeneous distribution and 
so their SF. 

Let us examine some details of the phase transitions 
in FigUl First, we focus on the transition from the PM 
phase to the AF phase that takes place as the value of ci 
is increased. In Figl^lwe present [/, C and p along C3 ~ 
4.0, 6 ~ 0.10 to locate the transition point between the 
PM and AF phases. It seems that the phase transition 
is of second order because there are no hysteresis in U 
and the peak of C develops systematically as L increases. 
We obtained similar behavior of U and C in the PM-AF 
phase transition for cases of other values of 6 shown in 
FigUl 

Secondly, we examine the transition between the PM 
phase and the FM-I-SF phase. In Fig|3]we present U,C 
along ci = 1.5 and fie = 10.0, which indicate a second- 
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FIG. 3. E, C and p along C3 = 4.0 and 5 ~ 0.10. A tran- 
sition between the PM and AF phases exists at ci ~ 2.9. 
Corresponding value of chemical potential pc = 10. 



order transition at C3 ~ 9.0 between the PM and FM-f SF 
phases, p increases as C3 increases crossing the transition 
point. 

Thirdly, we examine the transition between the AF 
and AF+SF phases at (5 ~ 0.15. In FiglS] we show U, 
C, p for C3 = 24.0 and pc = 16.0. There are hysteresis 
curves in these quantities, which indicate a first-order 
phase transition. Similar first-order transitions are found 
also for 5 = 0.05 and 5 = 0.15 as indicated by solid curves 
in FigU We determined the transition point from the 
hysteresis data by averaging the two values of p for fixed 
/ic, Ci, C3, pi on the upper hysteresis curve and p2 on the 
lower hysteresis curve, so that 



aipi + a2P2 — S, ai + a2 = 1, 



(3.6) 
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FIG. 4. U, C along ci = 1.5 and pc = 10.0. A second- 
order transition between the PM and FM+SF phases exists 
at C3 ~ 9.0. 
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FIG. 7. The specific heat for each term, Cb, Ch along ca = 
16.0 and ^c ~ 10.0. They exhibit the corresponding peaks(or 
cusp) for the transition between the FM+SF and AF+SF 
phases and the transition between the AF+SF and AF phases. 



FIG. 5. U, C, p vs ci for ca = 24.0 and fj.c = 16.0 (L = 12). 
A first-order transition between the AF and AF+SF phases 
exists at ci ~ 8.2. 



to obtain the given value oi 6. U and C are averaged 
with these weights ai,a2. 

Finally, we examine the successive transitions from the 
FM+SF phase to the AF+SF phase and then to the AF 
phase. In FiglS] we present U, C along c^ — 16.0 and 
IJ^c — 10.0. We see that there is only one peak in C 
at ci ~ 5.0 that corresponds to the transition between 
the FM+SF phase and the AF+SF phase, the peak for 
the transition from the AF+SF phase to the AF phase 
being missing. (For identification of these phases, see 
later discussion.) However the hole density p changes its 
behavior at Ci ~ 5.9, while C exhibits only very small 
anomalous behavior(dip) at that point. 

To clarify the possible transition at ci ~ 5.9, we mea- 
sured the "separate specific heats" for the A'^ term and 



the Ah term separately. They are defined as 



a = ^ {{AD - {A.r) , 



Gh = ^ {{AD (Ah)^) 



(3.7) 



In FiglT] we present Cg and Ch separately. Both Ch 
and Cs exhibit a very sharp peak at ci ~ 5.9. This kind 
of "cancellation" of specific heats of each term in the 
total specific heat C has been observed sometimes. This 
stems from the fact that {Ai,) is an increasing function 
of ci whereas (Ah) is a decreasing function of ci. At 
ci ~ 5.9, cancellation between these two terms occurs 
and {All) becomes a smooth function of C3. The study 
of correlation function in the next subsection shows a 
possible phase transition between AF+SF and AF phases 
at ci ~ 5.9. 



Correlation functions 




6.0 c, 6.4 



In the previous subsection, we showed the phase dia- 
grams for various hole density. In order to understand 
physical properties of each phase, let us study some cor- 
relation functions. Because the spin operator Sx of atom 
is given by Eq. (|2.6p as 



1 



BlaBr 



1 



{I - (f>l(f>x)zlazx, (3.1 



we define the corresponding normalized classical atomic 
spin vector. 



Sx = {'^ - 4>x4'x)ix, 



(3.9) 



FIG. 6. U, C, p along C3 = 16.0 and pc = 10.0. A second- 
order transition between the FM+SF and AF+SF phases ex- 
ists at ci ~ 5.0. There is a dip in C at ci ~ 5.9, which 
suggests another transition (See Fig[7]and the text). 



using the pure 0(3) spin vector H^ = z^azx of Eq. p.2l) . 
We also introduce the atomic density Ux corresponding 
to Eq. (12211), 

nx = \-(\)x(t>x- (3.10) 



Then we measure the following correlation functions, 
G,(r) = — y 



iN 



■x-\-rfi. ^x / 1 



Gsir) 



Gn{r) 



Gs{r) 



X,fJ, 

1 



6A^(1 - S) 
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3N{1-Sy 



/ ^ \-^x^rii.(7-'^xa / 1 

X.}.L.<7 

j" / ^ \P'X+r^'^x) •) 



X.fl 



3N{l-Sy 



■E(^^ 



i>x 



(3.11) 



From above, they are the spin correlation of spinous, 
atomic correlation, density correlation, and the spin cor- 
relation of atoms, respectively. Their prefactors are cho- 
sen so that G{r) is normalized as G(0) = 1. 

In Table 2, each phase is characterized by the 
(non)vanishing off-diagonal long-range orders (LRO) 
G(oo) and/or staggered magnetizations G(oo), 



G(oo) = lim G{r), 
G(oo) = lim {-YG{r), 



(3.12) 



measured by these correlation functions. 

To show these correlation functions, we select a typical 
point in the C3-C1 plane for each phase as in Table 3. In 
Fig|S]we present the four correlation functions of (13. lip 
at the four selected points in Table 3 for three typical 
values oi 6 = 0.15,0.20 and 0.30. From this result, we 
identified each phase as in FigHJ In the FM-I-SF phase 
and the AF-t-SF phase, the region 0.08 ^ 6 < 0.2 requires 
a difficult fine tuning of iic- So, to obtain the correlation 
data for 6 = 0.15 in these two phases, we used the super- 
position p.6p by using the data for 6 ~ 0.08 and 6 ~ 0.2. 

The results in FigIS] shows that the spin correlation 
Ge{r) and the atomic-spin correlation Gs{r) exhibit sim- 
ilar behavior in all phases. They have staggered magne- 
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Table 2. (Non) Existence of LRO in each phase. "SM" 
implies a nonvanishing staggered magnetization, G(oo). 
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FIG. 8. Four correlation functions p. lip at the four selected 
points in Table 3 for S — 0.15,0.20 and 0.30. They are used 
to identify the phases of Fig[l] according to Table 2. 
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Table 3. Data points of correlation functions for each 
phase in the C3-C1 plane. 



tizations in the AF and AF+SF phases. The atomic cor- 
relation Gsif) has nonvanishing LRO in the FM+SF 
and AF+SF phases, which indicates the SF of atoms. 
In particular, in the AF+SF phase, Gsir) of even-even 
(odd-odd) sites is sightly different from that of even-odd 
sites, reflecting the AF nature. The atomic density cor- 
relation Gn{r) shows smooth and homogeneous distri- 
bution of atoms and holes at all the examined cases of 
ci,C3,(5, indicating nonexistence of the PS phenomenon 




IV. 



RESULTS OF MC SIMULATIONS OF 
MODEL I 



FIG. 9. Four correlation functions (|3.11|) at the long distance 
limit G(I//2) vs ci for C3 = 16.0 and ^c = 10.0. These val- 
ues correspond to Figs l6l7l Gb{L/2) shows a behavior of a 
second-order phase transition at ci ~ 5.9 ~ 6.0. 



In this section, we shall study Model I of (I2.22p . Physi- 
cal meaning of the pref actor P^ = l — (j)x4>x in Ag of (|2.22p 
is obvious, i.e., if the site x is occupied by a holon, the AF 
coupling on the link {x,x ± fx) is suppressed. This bond- 
breaking effect of holons in the AF background generates 
an effective attractive force favoring holon pairs sitting 
at NN sites. Actually, a NN holon pair breaks eleven 
AF bonds while two holons separated by more than NN 
sites break twelve AF bonds [2J). This effective force 
can be an origin of the superconducting NN holon-pair 
condensation|l(]|| and/or of a PS of holon-rich regions and 
holon-free (holon-poor) regions. In this section, we shall 
study how the phase diagram is influenced by this AF 
bond-breaking effect by holons, in particular, if there ap- 
pears a phase-separated state. 



in the present model. This point will be discussed rather 
in detail in the following section. 

To study the possible phase transition between AF+SF 
and AF phases, we measured the correlation functions 
across the relevant region of ci. In FiglH] we plot their 
values at long-distance limit, G{L/2) for C3 = 16.0 and 
^c = 10.0. Because Gb{L/2) decreases continuously as 
ci increases and vanishes at ci ~ 5.9 ~ 6.0, we judge 
that there is a second-order phase transition from the 
AF-I-SF phase to the AF phase. Other three values 
G{L/2) remain finite as expected from Table 2. Simi- 
larly, FigUni shows G(L/2) for cg = 16.0 and ^ic = 16.0. 
Again Gb{L/2) decreases and vanishes but discontinu- 
ously with hysteresis. This supports our phase diagram, 
Figlll with the first-order AF+SF o AFtransition at 
ci ~ 8.2. 
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FIG. 10. Four correlation functions (|3.11[) at the long distance 
limit G{L/2) vs ci for C3 — 24.0 and /Xc = 16.0. These values 
correspond to FigO Gb{L/2) shows a behavior of a first- 
order phase transition at ci ~ 8.2 (dashed line), which is 
determined from the specific heat of FigO using (|3.6p . 



A. Phase structure 

As for Model II, we studied Model I by means of the 
MC simulations. In FigsfTTl we show the obtained phase 
diagram in the C3 — ci plane for fixed values of ^c- There 
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FIG. 11. Phase diagram of Model I in the C3 — ci 
plane for fixed values of the chemical potential, fic = 
-1-1.0, 0.0, —1.0, —2.0, —5.0. Hole density p is also shown for 
corresponding regions. The solid curves show first-order tran- 
sitions and the dashed curves show second-order transitions. 
The AF+SF phase is missing. 
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exist the AF, PM and FM+SF phases as in Model II with 
low-hole doping. In the phase diagram, we also indicate 
the hole density S. We found that hole density in the AF 
phase is rather small compared with Model II. This is due 
to the factor P^ in the ci-term of (I2.22p . As explained 
in the beginning of this section, a hole sitting on the 
site X makes P^ — and suppresses the AF couplings 
around x. So hole doping costs larger energy in Model I 
than in Model II. In other words, as the hole density is 
increased, the AF state becomes more unstable because 
the AF couplings are weakened by the prefactor P^+^Px- 
By the numerical studies of Model I, we have concluded 
that the coexisting phase AF+SF of AF order and SF 
does not exists in Model I. 

Order of each phase transition is indicated in FigHT] In 
particular, the transition between the AF and FM+SF is 
of first order. This result is physically expected as these 
two phases have different LRO's in contrast. We also 
found that the spin correlation G^, atomic correlation 
Gb, etc. exhibit a similar behavior to those in Model II. 



B. Phase separation and hole distribution 

The PS is a phenomenon of inhomogeneous distribu- 
tion of holes in which the lattice is separated into hole- 
rich regions and spin-rich (hole-free) regions. To study 
the possibility of PS, we introduced the following two 
quantities, 

^c = -^^^ ^ ((^x+M(/'a;+M - P){4>x<l>x ^ p)) 

^ ^^;p2^{(i>x+t,(l)x+t,(lyx(t>x~ p^) ■ (4.1) 
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FIG. 12. A and Ac of (gH) in Model I for C3 = 1.0 and 
/ic = 0.0 vs Cl. They have peaks around the phase transition 
point between the PM and AF phases, c\ ~ 3.6. 



In FiglT2| we first show A, Ac and p of Model I vs ci 
for C3 = 1.0 and jic = 0.0. A phase transition from the 
PM to AF phases takes place at ci ~ 3.6. Both A and 
Ac exhibit peaks around the transition point, though the 
magnitude of Ac itself is very small for all ci's. This be- 
havior around the transition point can be understood as 
a result of formation of AF domains of finite sizes in the 
PM phase on approaching to the AF phase boundary. In 
the AF phase, A has a fairly large value, which indicates 
that localization of holes takes place there. However, the 
small Ac means that the PS does not take place in both 
AF and PM phases, p decreases rapidly in the AF re- 
gion because holes break the AF bonds around them as 
explained before, costing higher energy in the AF phase 
than in the PM phase. 

The above results should be compared with results in 
Model II. See FigHH for Model II with cz = 5.0, pc = 
10.0, in which the PMoAF transition takes place and 



The quantity A measures the fluctuation of hole den- 
sity (i>x4>x around its average p. As explained in (|2.12p . 
a physical state at a specific site a; is a superposition 
of the one-holon state |0) {Bxa\^) — 0) and the no- 
holon states B|^|0). When the system is with homoge- 
neous distribution of holes like the holon-condensed state, 
4>x4'x — /C ~ and A '^ 0. On the other hand, when the 
system enters into an inhomogeneous hole-localized state, 
't'x't'x — p ^ —p for hole-free sites and 4>x4>x — p ^ 1 — p 
for sites occupied by holes, and A develops from zero. 

Similarly, Ac measures correlation of fluctuations of 
hole densities at the NN sites. When the deviation of hole 
densities from their average, (fixfj^x — P, at NN sites have 
simflar values (same signs), Ac develops a value of 0(1). 
A hole-rich region contributes amount of ^ (1 ~ p)^ to 
Ac and a hole-free region contributes ^ p~. If (j)x4'x takes 
finite but random values site by site, Ac ~ although 
A ~ 0(1)- In the phase-separated state, both A and Ac 
should take values of 0(1). 
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FIG. 13. A and Ac in Model II for C3 = 5.0, pc = 10.0 vs a. 
The transition between the PM and AF phases takes place at 
Cl ~ 2.9. 
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FIG. 14. A and Ac in Model I vs C3 for a = 1.0, fic = 0.0. 
The PM-f^FM+SF phase transition takes place at ca ~ 5.5. 



FIG. 16. A and Ac in Model I vs Ci for cg = 10.0, /ic = 0.0. 
The FM+SF <->-AF phase transition take a place at ci ~ 4.2. 



the hole density p takes almost the same value with that 
in the AF phase of Model I. However, A and Ac are 
almost constant while ci is varied across the phase tran- 
sition from the PM to AF phases. In Model II, the spin 
dynamics does not influence the hole dynamics for suffi- 
ciently small S. Obviously, the PS does not takes place 
in Model II as in Model I. 

Let us next see the results around the phase transition 
from the PM to FM-I-SF phases. In FigHH we show 
the results of Model I for a = 1.0, a*c = 0.0. The hole 
density 6 is fairly large, and both A and Ac are decreasing 
functions of C3. This result is expected, as holes in the SF 
phase are distributed homogeneously through the whole 
system. The corresponding results in the Model II for 
ci = 1.0 and p,c = 10.0 are given in FigfTSl Though 
the above parameters have different values from those 
in Figim the same PMoFM-l-SF phase transition takes 
place at C3 ~ 8.5. A exhibits similar behavior to that of 
the Model I, whereas Ac has a small peak at the phase 
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transition point. 

Finally, we show A and Ac around the FM+SFf^-AF 
phase transition. See FigHHfor Model I with C3 — 10.0, 
fic = 0.0. Both quantities exhibit hysteresis loops around 
the transition point ci ~ 4.2, because the phase transi- 
tion is of first order. In the AF phase, holes are strongly 
localized, whereas they have a homogeneous distribution 
in the FM-f SF phase. Corresponding numerical results 
for the Model II with C3 — 16.0, Hc = 10.0 are given in 
FigllTl There exist no anomalous behavior at the transi- 
tion FM-^SFoAF-hSF at a ~ 5.0, because the SF takes 
place in both phases. On the other hand, A develops in 
the AF phase as the SF disappears there. 

To close this section, we compare our results with other 
works. Our conclusion that there are no PS states is 
compatible with the results of Ref.fql for the isotropic 
2D model at low T's, although our system is 3D at finite 
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FIG. 15. A and Ac in Model II vs C3 for ci = 1.0, fj.c = 10.0. 
The PM^->FM-|-SF phase transition takes place at C3 ~ 8.5. 



FIG. 17. A and Ac in Model II vs a for C3 = 16.0, fj,c = 10.0. 
The FM-l-SFf^AF-fSF and AF-|-SF<->AF phase transitions 
take place at ci ~ 5.0 and ci ~ 5.9, respectively. See FigO 
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T's. As mentioned in Sect.l, Boninsegni and Prokof ev|8| 
studied the 2D bosonic t — J model with anisotropic spin 
couphng {Jx,y = ctJz^a < 1) with < J/t < 1 at low 
T's, and found PS for low hole concentrations, S < 6c{a), 
where Sc{ct) is the critical density of holes and dc{a) — > 
as a — > 1. This implies that anisotropy is essential for 
PS; the stability of spins in z-direction via the NN spin 
coupling favors to form spin-rich regions. 



V. CONCLUSION AND DISCUSSION 

In the present paper, we proposed the bosonic t-J 
model of hard-core atoms and studied its physical prop- 
erties by means of the MC simulations. In Model II, 
the effects of the projection operator Px — 1 — (j)l.(j)x are 
ignored and, as a result, there appear four phases, AF, 
PM, PM-hSF, and AF-f SF phases. In Model I, we mostly 
studied the bond-breaking effects by holes and found that 
the AF-I-SF does not appear in the phase diagram. We 
also investigated if the PS takes place and found that, in 
both models, it does not. 

We argued that Model II describes effectively a doped 
AF Heisenberg model including strong AF correlations 
besides the NN AF couplings, such as the next-NN spin 
couplings. Thus the important result of Model II, appear- 
ance of the AF-t-SF phase, suggests that such a coexisting 
phase of AF order and SF may be generated by inclusion 
of certain AF-enhancing couplings such as the next-NN 
coupling beyond the NN ones. 

For the fermionic t-J model, whether the AF+SF 
phase (or rather the coexisting phase of AF order and su- 
perconductivity) exists or not is an open question. How- 
ever a AF-|-superconducting phase has been observed by 
recent experiments with uniform high-quality samples of 
high-T cuprates[25|. Therefore we expect that a similar 
coexisting phase of AF-f SF appears in cold-atom sys- 
tems in an optical lattice with strong on-site repulsion as 
long as certain AF-enhancing couplings such as next-NN 
coupling are involved. 

Another interesting subject is to consider the effect of 
anisotropy of pseudo-spin coupling {a < 1). Although we 
studied the isotropic case (a = 1) in the present paper as 
the canonical case, the realistic cold-atom systems may 
have anisotropy. This is because the difference in two 
species of bosons such as their hopping amplitudes, etc., 
are reflected to this anisotropy in the resulting t-J model 
(recall the relation J = 4t'^/U). 

At the beginning of Sect. 4, we mentioned the possible 
superconductivity via formation and condensation of NN 
hole pairs in Model I. Because we found no separated 
phase transition into superconducting phase in Sect. 4, it 
is plausible that the superconducting state of hole pairs 
is generated simultaneously with the SF state of single- 
boson EEC. To examine this point further, it is necessary 
to study the correlation function of these pairs in de- 
tails. When the anisotropy is included, such correlation 
may have different behaviors in different pairing chan- 



nels. These topics (the next-NN coupling, anisotropy, 
and the hole-pairing) are under study. We hope to re- 
port on them in future publications. 

Finally, studies on the bosonic t-J model in the 
present and previous papers suggest for the fermionic 
t-J model that, at finite hole concentrations, coherent 
holon hopping is realized. In the slave-fermion approach, 
this implies that a small Fermi surface of fermionic 
holons is generated in under-doped region. In the 
further study of the bosonic t-J model, we expect to 
draw further interesting suggestions for the fermionic 
t-J model and vice versa. 
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Appendix A: CP^ Representation of a Hard-core 
Boson 

In this Appendix we derive the representation (jl.3[ 
12.181) of a HCB operator in terms of a CP^ spin oper- 
ator Wxri- 

Let us introduce an s = 1/2 SU(2) (pseudo) spin op- 
erator Tx by using Wxri as 

Tx = -Wx^Wx, (A.l) 

where a are 2x2 Pauli matrices. Explicitly, 

Txl = -^{WxlXx2 + Wx2Xxl), 

Tx2 = --^iWxlXx2 - Wx2Xxl), 

Tx3 = -^{WxlXxl - Wx2Xx2)- (A. 2) 

On the other hand, the rising operator T^" and the 
lowering operator T,^ defined by 

(A.3) 

(A.4) 



T^ = Tix ± iT2x, 



is expressed by the HCB (px as 
Then we have 



Tx = rx. T^ = <l>x- 



<l>x = 'WI2WX1, (A. 5) 

which is (|1.3[|2.1g)) . The correspondence between the two 
sets of states are as follows; 

^31+)-^!+), Ts\-) = -^\-), 

0|O)=O, |-) = |0), |+)=0^|O), 
Wn\v)=0, |-)=w^|i;), \+)^wl\v). (A. 6) 
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